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INTRODUCTION 


During /iperiod  13  September  1985  to  31  March  1987, 

KM  Sciences  -furnished  support  to  the  radiation  e-ffects 
programs  of  the  Radiation-Matter  Interactions  Branch  of  the 
Naval  Research  Laboratory  in  the  areas  of  collecting, 
manipulating,  analyzing,  and  displaying  experimental  data, 
in  writing,  modifying,  converting,  and  extending  computer 
codes  for  modeling  physical  phenomena  and  in  graphical 
presentation  of  experimental  data  and  results  of 
computations.  Two  major  tasks  were  supported: 

(1)  computer  simulation  of  the  effects  of  particle  beams  on 
materials  and  (2)  experimental  measLu^ements  of  the  effects 
of  particle  beams  on  materi^s.  '  NffdPoP  Cfy4yJs>^ 

ACCOMPLISHMENTS 


COMPUTER  SIMULATION: 


MARPOP,  a  computer  program  (written  for  the  TI-ASC 
computer)  that  processes  output  from  the  MARLOW  cascade 
simulation  code,  was  converted  to  run  on  NRL's  Cray 
computer.  A  new  output  routine  for  the  Cray  version  of 
MARLOW  was  written  to  produce  an  output  data  file  that 
serves  as  input  to  MARPOP.  MARPOP  performs  calculations 
and  manipulations  on  the  data  from  MARLOW  and  produces 
tables  and  graphs  summarizing  and  displaying  the  MARLOW 
results.  A  number  of  production  runs  were  made  with  the 
Cray  versions  of  MARLOW  and  MARPOP  using  input  data 
designed  to  simulate  experiments  in  which  materials  were 
irradiated  with  particle  beams.  In  an  evolutionary 
process,  the  experimental  data  and  the  results  of  the 
computer  simulations  led  to  frequent  changes  in  the  MARPOP 
processing  and  output  specifications,  requiring  numerous 
changes  in  the  MARPOP  code.  The  results  of  these 
calculations  were  incorporated  in  a  paper,  "Cascade 
Simulation  of  the  Crystal  Orientation  Dependence  of 
Sputtering  and  Lattice  Damage  of  Single  Crystal  Copper  by 
Irradiation  with  100  keV  Copper  Ions",  by  8. P. Mueller, 

M. Rosen,  W. A. Fraser  (KM  Sciences),  J. A. Sprague, 

P. R. Malmberg,  J.M. Lambert,  P.A.Treado,  and  G. W. Reynolds, 
published  in  Nuclear  Instruments  and  Methods  in  Physics 
Research  818(1987)360-364.  A  source  code  listing  of  the 
most  recent  version  of  MARPOP  is  enclosed. 


PARTICLE  BEAM  EXPERIMENTS: 


Computer  support  was  provided  for  five  series  of 
experiments  using  the  NRL  Linear  Accelerator  to  study  the 
effects  of  particle  beams  on  materials.  The  support 


1 


con«ist«<j  of  computor  data  acquisition,  data  manipulation 
and  calculations,  and  graphical  and  tabular  display  of 
rosults. 

OTHER: 

A  hands-on  guide  with  sample  command  sequences  was 
developed  to  introduce  users  to  graphics  software  installed 
on  the  Condensed  Matter  and  Radiation  Sciences  Division  VAX 
computer.  This  enabled  members  of  the  Radiation-Matter 
Interactions  Branch  to  produce  many  types  of  data  plots 
with  less  than  an  hour  of  self-training  using  the  guide. 

Hardware  and  software  specifications  were  developed 
for  a  personal -computer-based  pilot  system  for  storage  and 
retrieval  of  experimental  data  and  results  of  theoretical 
calculations  of  effects  of  charged  particle  beams  and 
deposition  of  energy  in  materials.  After  the  personal 
computer  systems  were  delivered,  they  were  set  up  and 
software  was  installed. 
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PROGRAM  MARPOP 


VERSION  3.  1  22  MAY  1986 

Programnvr:  U.  A.  FRASER  KM  Scisncv* 

MARPOP  (MARlow  P0»t-PTOC»5»or )  reads  Cray  program  rVkRLOW  'POPDAT' 
output  fils  (data  for  each  particle  exiting  from  the  target  for  zero 
surface  binding  energy)  and  recalculates  the  yield*  energy  and 
directional  distributions  for  chosen  surface  binding  energies. 

THIS  VERSION  OF  MARPOP  PROCESSES  OUTPUT  FILES  FROM  MARLOW  VERSION  OF 
18  (+  OR  -1)  FEBRUARY  1982*  ADDING  TABLES  OF  MAXIMUM  PARTICLE  DEPTH* 
AND  NUMBERS  AND  LENGTHS  OF  COLLISION  SEQUENCES  TO  THE  OUTPUT  THAT  WAS 
PRODUCED  BY  MARPOP  VERSIONS  1.3  AND  1.4.  THE  "18**  FEBRUARY  MARLOW 
VERSION  TRUNCATES  MAXIMUM  DEPTH  VALUES  TO  THE  NEXT  LOWER  0.  1  ALATd) 
UNIT.  DATA  FROM  UP  TO  10  MARLOW  OUTPUT  FILES  CAN  BE  COMBINED.  THE 
NUMBER  OF  PRIMARY  PARTICLES  FOR  WHICH  DATA  WILL  BE  PROCESSED  FROM  A 
FILE  CAN  BE  SET  TO  LESS  THAN  MAXRUN.  SURFACE  BINDING  ENERGY  FOR 
PRIMARY  PARTICLES  IS  0.  0. 

*****  Input  Parameter  Records  ***** 

*****  read  from  logical  unit  LUSPEC  ***** 

Up  to  20  parameter  sets  may  be  input.  Each  set  consists  of  the 
following : 

Record  Format 

1  Title  may  contain  any  readable  characters  A80 

(maximum  *  80  characters) 

2  hFlLES  lOX.  12 

Number  of  MARLOW  ^POPDAT'  files  to  be  processed 
(maximum  >  10). 

3  LIMRUN(NFILES)  lOX*  1015 

f^FILES  values  of  LIMRUN*  the  number  of  primary 
particles  for  which  data  are  to  be  processed  from 

each  'POPDAT'  file.  The  first  value  on  record  3 
will  be  used  for  file  'FT17'*  the  next  for  'FT18'* 
etc.  If  LIMRUN(i)»0*  the  corresponding  FTnn  file 
will  be  processed  to  the  end  (MAXRUN  primaries). 

If  LIMRUN(i)>0«  then  data  from  the  corresponding 
FTnn  file  will  be  processed  until  data  produced  by 
primary  particle  LIMRUN4>1  are  encountered. 

4*5  SBND(l-lO)  lOX,  5E10  0 

Surface  binding  energies  (e.  v.  )  for  up  to  10 
particle  types.  MARLOWE  Version  12  provides  for 
only  5  particle  types.  Ten  types  were  permitted 
in  Version  11.  Both  records  4  and  5  must  be 
present  (to  process  MARLOW  Version  12  output  there 
will  be  no  values  on  record  5). 

6*7  WIDTH(l-lO)  lOX,  5E10  0 

Widths  (e.  V.  )  of  channel*  for  binning  the  energies 
of  up  to  10  types  of  secondary  particles.  MARLOWE 
Version  12  provide*  for  only  5  particle  types  Ten 
types  were  permitted  in  Version  11.  Both  record* 


OOOOOOOOOOOOOOOOOOOOOOOOOOOrtOrtOOrtOOOOOnrtriOOOOOOOOOOOOOOOOOOO 


4  Afitf  7  Mi«t  b«  pr*s*nt  (to  procooo  HMIL(X4  Voroion 
12  output  thoro  will  bo  no  voluot  on  rocord  7). 

t  yXDTHdl)  lOX.  EiO.  0 

Width  (o.  V.  )  of  chonnols  for  binning  tho  onorgio* 
of  priMorg  port  id  os. 

9  DEPBZN.  DEFHAX  lOX,  2E10.  0 

DEPBIN  >  bin  width  (in  units  of  ALAT(l))  for 
dopth  distribution  tab  Its. 

OEPMAX  «  bin  width  (in  units  of  ALAT(1>)  for 
moiiisuiB  dopth  distribution  toblos. 

10  NnUBXN.  CHnU  lOX.  12.  1X.A2 

NMUBIN  nuMbor  of  mu  (cosino  polar  anglo)  to  uso 
for  binning  particlos.  Prosont  maximum 
is  20.  limitod  by  procossing  and  output 
routinos.  but  input  routino  providos  for 
40  bins.  Arrays  in  othor  routinos  and 
output  formats  must  bo  modifiod  bofore 
40  bins  can  bo  usod. 

CHMU  «  data  typo  spocifior.  OP1U«'nU'  if  bin 
boundary  valuos  on  rocord(s)  11  aro  in 
units  of  mu  (cosino  polar  anglo).  Loavo 
CHMU  blank  if  boundary  valuos  aro  in 

dogroos. 

IIA.  IIB.  otc.  aMU(0  to  NMUBIN)  or  P0Ll>Ee(0  to  NMUBIN)  lOX.  7F10.  0 

As  many  rocords  as  aro  nocossary  for 
NMUBIbH'l  valuos  of  bin  boundarios  for  binning 
polar  anglo  (thota)  distribution  of  particlos. 

If  CHMU  (rocord  10)  is  bin  boundaries  must 

bo  cosinos.  If  CHMU  is  anything  olso.  bin 
boundarios  must  bo  in  dogroos. 

12  NAZBIN  lOX. 12 

Number  of  azimuthal  anglo  (phi)  bins  to  uso  for 
binning  particlos.  Prosont  maximum  is  20. 
limittod  by  procossing  and  output  routinos.  but 
input  routino  providos  for  40  bins.  Arrays  in 
othor  routinos  and  output  formats  must  bo  modifiod 
boforo  40  bins  can  bo  usod. 

13A.  13B.  otc.  AZK)E0(0  to  NAZBIN)  lOX.  7F10.  0 

As  many  rocords  as  aro  nocossary  for  NAZBIN-*-! 
valuos  of  bin  boundarios  for  binning  azimuthal 
anglo  (phi)  distribution  of  particlos.  Values 
must  bo  in  dogroos. 

14  XN0RI9.  lOX.  FIO.  0 

Factor  (mmjI  tipi  ior)  to  bo  usod  in  normalizing 
valuos  for  ojoctod  particle  direction  distribution 
table.  Boo  writeup  describing  normalization  of 
array  NANOLE  in  subroutine  SUMRYZ. 

19  OUTPUT  CONTROL  SWITCHES  lOX. A70 

To  eliminate  undosirod  output,  one  or  more  of  tho 
following  codes  (separated  by  commas  or  spaces) 
may  appear  in  columns  11-80.  in  any  order. 

FPRIM  ("Front  PRIMarios")  supprossos  output  for 


ooofto  o  n  oo  n  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  rt  oorto 


r«#l«ct»4  prinarids. 

RPRZH  ("Rdar  PRXHariPs")  tupprppsdft  output  for 
tranoaittod  priaariao. 

fn'AflO  ("Front  TAROat")  auppraaaaa  output  for 
front  aputtarad  tar  gat  atoaa. 

RTARO  ("Raar  TAR Oat " >  auppraaaaa  output  for 
raar  aputtarad  targat  atoaa. 

XfFO  cauaaa  fila  'XNFO'  not  to  ba  includad  at 

tha  and  of  tha  print  fila-  Tha  INFO  fila 
aag  ba  aavad  or  diapoaad  of  juat  aa  any 
*othar  fila. 

Tba  racord  auat  ba  praaant  ao  if  all  output  ia 
daairad«  coluana  11-80  ahould  ba  blank. 

aaaaa  OTHER  INPUT  a#*** 


Data  filaa  (naaad  TOPOAT')  writtan  by  prograa  HARLOW  froa  routina 
EXTRAX.  Froa  1  to  10  POPOAT  filaa  aay  ba  procaaaad  in  a  aingla 
HARPOP  run.  All  POPOAT  filaa  ara  procaaaad  for  aach  HARPOP  input 
paraaatar  aat  daacribad  abova.  POPOAT  filaa  ara  raad  froa  logical 
unita  17i  18* .  .  . .  «  25<  2A«  and  auat  ba  aaaignad  tha  naaaa  FT17*  FT18> 
ate.  I  up  to  FT26. 

ETC.  ). 


##»*«  VARIABLE  DECLARATIONS  aaaaa 

Carriage  control  charactara 
CHARACTER* 1  CHPLUS.  CHWUN 
Output  awitchaa 

CHARACTERaS  FPRINi  RPR  IN.  FTARO.  RTARO.  INFO 
CHARACTERaS  CHKQDE 
CHARACTERaeO  CHKARD 
CHARACTER* 130  CHLINE 

INTEGER  X.  18ET.  KODE.  LOATE.  LIHRUN.  LTIHE.  LUNIT. 

1  HXFYLS.  NAZBIN.  fFILES 

REAL  DEPBIN.  DEPHAX.  SBND.  WIDTH 

PRIHARY  PARTICLE  PROCESSING  LIMITS  READ  FROM  CARD(S>  3 
COMMON  /LIMITS/  LIMRUN(IO) 

Output  awitchaa 

COMMON  /OUTSWT/  FPRIM,  RPRIM,  FTARO.  RTARO.  INFO 
COMMON  /READIN/  SBND(ll).  UIDTHdl).  DEPBIN,  DEPMAX 
I/O  logical  unita  aat  in  aubprogram  BLOKAA 
COMMON  /UNITS/  LUSPEC.  LUPRT.  LUINFO 
NEXT  STATEMENT  FORCES  LINKING  OF  BLOCK  DATA  SUBPROGRAM  THAT  SETS 
LOGICAL  UNIT  NUMBERS. 

Naxt  atataaant  forcaa  linking  of  block  data  aubprogram  that  sata 
conatant  valuaa  for  I/O  logical  unita.  ate. 

EXTERNAL  BLOKAA 

DATA  CHPLUS  /'♦'/,  CHWUN  /  'I'/ 

C  Maximum  numbar  of  POPDAT  filaa  that  may  ba  procaaaad 

DATA  MXFYLS  /lO/ 

C 

100  FORMAT (ABO) 

101  FORWkTdOX.  5E10.  0) 

102  FORMAT  ( lOX.  2E 10.  0.  IS) 

103  FORMAT (lOX. 12) 


104  FORMAT (lOX I  1013) 

105  FORMAT (A130) 

200  FORMAT (Al/ 'O'.  7X.  'PROGRAM  MARPOF  -  MARLOW  POST  PROCESSOR '.  15X. 
1  'VERSION  3.  1'.  5X.  'OF  22  MAY  1986') 

201  FORMAT<'  '.TX.  '***  UNEXPECTED  END-OF-FILE.  INPUT  DATA  CARDS') 

202  FORMAT ('  '.7X.  '»**  READ  ERROR.  INPUT  DATA  CARDS') 

203  FORMAT <'  '.7X.  'PROCESSING 14.  2X,  'POPDAT  FILES') 

204  FORMAT<'  '.7X.  'PROGRAM  MARPOP  FINIS«D') 

205  F0RMAT</'0'.27X.  A80/) 

206  FORMATC'  '.7X.  '***  UNEXPECTED  END-OF-FILE.  UNIT  '.  14) 

207  FORMAT('  '.7X,  '***  READ  OR  DATA  ERROR.  UNIT  '.  14) 

208  FORMAT  ('  '.7X. 

1  '«««  MAXIMUM  NUMBER  OF  RECORDS  PROCESSED  FROM  UNIT  '.  14) 

209  FORMAT ('  '.7X,  '***  END-OF-FILE  READ.  INPUT  DATA  CARDS') 

210  FORMAT  ('  '.7X.  'UNIT  '.  14.  3X,  'MAXRUN  PRIMARIES  TO  BE  PROCESSED') 

211  FORMAT<'  '.7X.  'UNIT  '.  14,  1 10,  2X.  'PRIMARIES  TO  BE  PROCESSED') 

212  FORMAT  < 'O'.  7X,  'DATE ',  1  X.  A8.  3X.  'TIME '.  1  X.  A8  ) 

213  FORMAT { 'ICONTENTS  OF  FILE  INFO') 

214  FORMAT (A1 30) 

215  FORMAT('  LAST  RECORD  OF  FILE  INFO') 


C 

C 

C 

C 

C 

C 


C 


C 

C 

C 

c 

c 

c 

c 

c 


WRITE  (LUPRT.  200)  CKPLUS 
CALL  DATE  (LDATE) 

CALL  CLOCK  (LTIME) 

OPEN  (UNIT-LUINFO.  FILE* 'INFO'.  STATUS* 'NEW' > 

WRITE  (LUPRT.  212)  LDATE.  LTIME 
WRITE  (LUINFO.  200)  CHPLUS 
WRITE  (LUINFO.  212)  LDATE.  LTIME 

On  »mch  past  through  DO  4000  ISET  loop,  one  set  of 
parameter  records  is  read  from  logical  unit  LUSPEC  and 
proc  essed. 


DO  4000  ISET-1.20 

*♦*  Record  1  *** 

READ  (LUSPEC,  100.  END*8000.ERR*10000>  CHKARD 
IF  (ISET.  GT.  1)  WRITE  (LUPRT,  200)  CHWUN 
WRITE  (LUPRT,  205)  CHKARD 
WRITE  (LUINFO.  205)  CHKARD 


Record  2  **•* 

READ  (LUSPEC.  103.  END*9000,ERR»1 0000)  NFILES 
IF  (NFILES.  GT.  MXFYLS)  NFILES-MXFYLS 
WRITE  (LUPRT.  203)  NFILES 

***  Record  3  *** 

READ  (LUSPEC.  104,  END-9000,  ERR-10000)  (LIfRUN(  I ).  I-l .  NFILES) 


♦**  Records  4,5  *** 

READ  (LUSPEC.  101.  END-9000.  ERR-10000)  (SBND(  I ),  I-l.  10) 


**♦  Records  6.7.8  *** 

READ  (LUSPEC.  101.  END-9000.  ERR-10000)  WIDTH 

**♦  Record  9  *** 

READ  (LUSPEC. 102. END-9000.  ERR-10000)  DEPBIN.  DEPMAX 


Remaining  parameter  records  are  read  from  LUSPEC  in 
subroutine  INPTB. 

CALL  INPTB  (CHKODE) 

IF  (CHKODE.  EQ.  'EOF  ')  THEN 

GO  TO  9000 

ELSE  IF  (CHKODE.  EQ.  'ERROR  ')  THEN  I 

GO  TO  10000  I 


ooooooooo  oooooo  ooo  f>ooo 


END  IF 
CALL  INZTLZ 

On  ffach  pass  through  00  3000  JFYL  loop  ons  MARLOW 
'POPOAT'  fils  is  road  and  particls  data  tabulatod. 

DO  3000  JFYL-i.  NFILE8 
LUDAT  -  16  +  JFYL 
IF  (ISET.  OT.  1)  REWIND  LUDAT 
IF  (LIMRUN(JFYL>.  LE.  0)  THEN 
WRITE  (LUPRT.  210)  LUDAT 
gLBE 

WRITE  (LUPRT.  211)  LUDAT.  LIMRUN(JFYL) 

END  IF 

Road  first  record  from  logical  unit  LUDAT 
CALL  INPTA  (JFYL.  LUDAT.  KODE) 

Values  of  KODE  returned  by  INPTA; 

1  a  Data  road  OK 

2  "  End-of-file  read,  no  FTxx  file 

3  a  Read  error,  file  FTxx 

GO  TO  (1000.  9000.  7000)  KODE 
1000  CONTINUE 

CALL  PRTCLS  (ISET.  JFYL.  LUDAT.  KODE) 

VALUES  OF  KODE  RETURNED  BY  PRTCLS: 

1  RECORDS  READ  AND  PROCESSED  OK 

2  UbEXPECTED  END-OF-FILE.  NO  DATA 

3  UNEXPECTED  END-OF-FILE.  PARTIAL 

DATA 

4  READ  ERROR.  NO  DATA 

9  READ  ERROR.  PARTIAL  DATA 

6  MAXinun  NUMBER  OF  RECORDS  READ 

AND  PROCESSED 

GO  TO  (30(X}.  9000.  6000.  7000.  10900.  2000)  KODE 
2000  CONTINUE 

WRITE  (LUPRT. 208)  LUDAT 
3000  CONTINUE 
CALL  SUMRYZ 

C  PRINT  TABLES  AND  HISTOGRAMS 

CALL  OUPTAA 
4000  CONTINUE 

GO  TO  12000 
9000  CONTINUE 

WRITE  (LUPRT.  206)  LUDAT 
IF  (I.  OT.  I )  GO  TO  11000 
GO  TO  12000 
6000  CONTINUE 

WRITE  (LUPRT.  206)  LUDAT 
GO  TO  11000 
7000  CONTINUE 

WRITE  (LUPRT.  207)  LUDAT 
IF  (I.  OT.  1  )  GO  TO  11000 
GO  TO  12000 
8000  CONTINUE 

IF  (ISET.  LE.  1)  GO  TO  9000 
WRITE  (LUPRT. 209) 

GO  TO  12000 


9000  CONTINUE 

WRITE  (LUPRT.201) 

00  TO  12000 
.10000  CONTINUE 

WRITE  (LUPRT,202) 

GO  TO  12000 
10900  CONTINUE 

WRITE  (LUPRT<207}  LUDAT 
11000  CONTINUE 

C  PRINT  TABLES  AND  HISTOGRAMS 

CALL  OUPTAA 
12000  CONTINUE 

WRITE  (LUINFO.  219) 

ENOFILE  LUINFO 
WRITE  (LUPRT.  204) 

IF  (INFO.  EQ.  'INFO  ')  THEN 
REWIND  LUINFO 
WRITE  (LUPRT.  213) 

DO  13000  1*1.600 

READ  (LUINFO. 109. END*14000)  CHLINE 
WRITE  (LUPRT. 214)  CHLINE 
13000  CONTINUE 
END  IF 

14000  CONTINUE 
C 

C  Next  stateoient  is  only  normal  program  execution  stop 

C 

STOP 

END 

SUBROUTINE  ANOPLT  (KTYPE.  KSURF.  AVYELD.  ERROR) 

C 

C  13  MARCH  1986 
C 

C  ANOPLT  MAKES  LINE  PRINTER  HISTOGRAM  PLOT  OF  ANGULAR  DISTRIBUTION 
C 

C  »*«*«  CALLING  PARAMETERS  *«*«« 

C 

C  KTYPE  -  (INTEGER.  PASSED)  PARTICLE  TYPE  NUMBER  (1-10  FOR  SECONDARY, 

C  11  FOR  PRIMARY  PARTICLES) 

C  KSURF  -  (INTEGER.  PASSED)  TARGET  SURFACE  ID  (1-FRONT.  2-BACK) 

C  AVYELD  -  (REAL.  PASSED)  HEM4  SECONDARY  PARTICLE  YIELD  PER  PRIMARY 
C  ERROR  -  (REAL.  PASSED)  STANDARD  DEVIATION  OF  fCAN  YIELD  (AVYELD) 

C 

C  THE  PLOT  IS  SET  UP  IN  <WRAY  LYhPLT,  81  COLUMNS  (1-LEFTMOST)  BY  42 
C  LINES  (LINE  2  IS  TOP  OF  PLOT  FRAME.  LINE  1  IS  PRINTED  ABOVE  THE 
C  FRAME).  IF  COSINE  BIN  20  (NPLOT( 20.  KTYPE.  KSURF) )  IS  GREATER  THAN 
C  ZERO.  THE  HEIGHT  OF  THE  RIGHTMOST  HISTOGRAM  BAR  -  1.0.  AND  THE  OTHER 
C  BAR  HEIGHTS  ARE  PROPORTIONED  TO  IT.  IF  NPLOT (20.  KTYPE.  KSURF )  -  O, 

C  THE  LARGEST  BIN  IS  FOUND  AND  ASSIGNED  A  HEIGHT  OF  1.0,  AND  THE  OTHER 
C  BARS  ARE  PROPORTIONED  TO  THAT  ONE.  IF  A  BAR  HAS  A  CALCULATED  HEIGHT 

C  GREATER  THAN  1.0.  THE  TOP  OF  1TE  BAR  EXTENDS  JUST  ABOVE  THE  PLOT  AND 
C  THE  CALCULATED  HEIGHT  (ENCODED  INTO  ARRAY  IXCEED)  IS  PRINTED  ABOVE 
C  THE  TOP  OF  THE  BAR. 

C 

C  ***♦♦  VARIABLE  DECLARATIONS  ***** 

C 

C 

C  PLOT  AXIS  AND  HEADING  LABELS 

DIMENSION  LABLAX(ll).  LAB8RF(4. 2) 

C  STORAGE  FOR  PLOT  CHARACTERS 


OO  O  ft  o  o  o  ooooo  o  oo 


c 

c 


DXHENBZQN  1XCEED(20).  LYNPLT(81.  42).  LYN0UT(324.  10 ) .  LE01(24). 

1  LE02(24).  LE03(24).  LE04(24) 

TEMPORARY  STORAGE  FOR  PLOT  LEGEND  LINES 
DIMENSION  LEOENDO) 

COMMON  /BINS/  P0L0E0(0: 40) .  QMU(0:40).  DMU(40).  NMUBZN. 

1  AZMDE0(0:40).  AZMRAD(0: 40).  0AZBIN<40).  NAZBIN 

ADUrtlY  «  ARRAY  NOT  USED  IN  THIS  SUBPROGRAM 
LTYPE  »  PARTICLE  ID  (LITERAL) 

COMMON  /FSTREC/  A0UmY(4a).  LTYPE(IO).  DDUMMY(70) 

SURFACE  BINDING  ENERGIES 
COMMON  /READZN/  SBNO(ll) 

ENERGY  (NENERG).  YIELD  (NYIELO).  DEPTH  (NDEPTH).  AND 
ANGULAR  (NANGLE)  DISTRIBUTION  TABLES. 

KSUM.  KSUMSQ  -  SUM  AND  SUM  OF  SQUARES  OF  YIELD  FOR 

CALCULATING  MEAN  YIELD  AND  STANDARD  DEVIATION. 
NPRYMS  >  TOTAL  NUMBER  OF  PRIMARY  PARTICLES. 

COMMON  /TABLES/  NENERG  ( 100.  11.  2).  NYIELD(21.  10.  2) . 

1  NDEPTH (21.  10.2).  NRn.ZD(21. 21.  11.  2).  FACNRM(  1 1 . 2) . 

2  NANGLE (21.  21.  11.2).  KSUM(10.2).  KSUMSQ ( 10.  2). 

3  NPRYMS 

COMMON  /UNITS/  LUSPEC.  LUPRT.  LUINFO 


B 

A 

M 

1 

m 

I 


EQUIVALENCE  (LEGl  ( 1 ).  LYhPLT(4.  4)  ).  (LE02(  1  ) .  LYNPLT  (4.  9  )  ). 

1  (LE03(  1).LYNPLT(4.  6)).  (LE04(1  ).  LYNPLT(4.  7) ). 

2  (LYNOUTd.  D.LYWLTd.  2)) 

SAVE 


DATA  LABLAX 

/'O.O  '. 

'0.  1  '. 

'0.2  '. 

'0.3  ', 

'0.4  '. 

'0.5  ', 

1 

'0.6  '. 

'0.7  '. 

'0.8  '. 

'0.9  '. 

'1.0  '/ 

DATA  LABSRF 

/'  RE', 

'FLEC  ', 

'  FR '. 

'ONT  '. 

'TRAN'. 

'SMIT'. 

1  '  R',  'EAR  '/ 

NUMBER  OF  CHARACTERS  IN  ARRAY  LEGEND 
DATA  NLEG  /28/ 

SYMBOLS  FOR  PLOT 

DATA  I  BLANK  /'  '/.  IDASH  /'-V,  IPLUS  /'*'/ 

DATA  MXAZM  /21/ 

201  FORMAT  ( IHl //// IHO,  34X.  'ANGULAR  DISTRIBUTION'.  3X,  IX.  2A4. 

1  'SPUTTERED  '.  2X.  A2.  13.  2X,  'ATOMS  ' ) 

202  FORMAT (IHl//// IHO,  33X.  'ANGULAR  DISTRIBUTION '.  3X.  3X.  2A4. 

1  'TED  PRIMARY  PARTICLES') 

203  FORMAT  (IH  .24X.20A4) 

204  F0RMAT((1H  .  19X.  A4.  81A1.  IX.  A3/  2dH  .23X.  81Al/>.  IH  .23X.  81A1)) 

205  FORMAT (20X.  21  (3X.  'I')  /  22X,21(F4.  1)  / 

1  '0'.49X.  'Polar  Anglo  (Dogroos)'/) 

206  FORMAT < IHO. 7X.  '***  NO  DATA  FOR  ANGULAR  DISTRIBUTION  PLOT'/  IHl) 

207  FORMAT(  'O'.  19X.  '1.  0-',  18.  2X.  'Part i c loo/dMu  ' ) 

401  FORMAT  (F4.  1) 

402  FORMAT  (  'SURF.  BINDG.  EN.  '.  1PE9.  2.  IX) 

403  FORMAT  ( 'PRIMARIES'.  9X.  17,  3X) 

404  FORMAT  ('MEAN  Y  lELD  ' ,  IX .  FIO.  1,  3X) 

405  FORMAT  ( 'ERROR',  6X,  FIO.  1.3X) 


m 


g 


c 


IF  (KTYPE.LT.il)  PRINT  201. 
1  KTYPE 

IF  (KTYPE.EQ.il)  PRINT  202. 
00  800  1-1.20 


PRINT  HEADING 

(LABSRFd,  KSURF).  1-3.  4).  LTYPE  (KTYPE) , 

(LABSRFd.  KSURF).  I-l.  2) 

BLANK  PLOT  ARRAY.  PUT  IN  GRID  LINES 


IXCEED(I)  •  IBLANK 
800  CONTINUE 

DO  900  1-1.42 
DO  900  J-1.  81 

LYNPLT(J. I)  -  IBLANK 
900  CONTINUE 

DO  1000  1-2.42.8 
DO  1000  J-1. 81.  2 

LYNPLT<J. I)  -  IDASH 
1000  CONTINUE 

DO  2000  J-2.  42 

DO  2000  I-l.  81.  16 

LYNPLT(I.J)  -  IDASH 
2000  CONTINUE 

C  SET  HISTOGRAM  BAR  HEIGHT  REFERENCE 

lYTOP  -  NRMLZDd.MXAZM.KTYPE.  KSURF) 

DO  3000  I-NMUBIN.  NtIUBlN 

IF  (NRMLZDd.  MXAZM,  KTYPE.  KSURF).  GT.  I  YTOP  ) 

1  I  YTOP»NRMLZD<  I.  MXAZM.  KTYPE,  KSURF) 

3000  CONTINUE 

IF  (lYTOP.  GT.  0)  GO  TO  4000 
PRINT  206 
RETURN 


4000  CONTINUE 

YMAX  -  FLO AT < I YTOP) 

LYNPLTd.42)  =  IPLUS 
last  »  42 


BET  UP  HISTOGRAM  IN  PLOT  ARRAY 


5000 


6000 


7000 

8000 

9000 


HISTOGRAM  BAR  FDR  ONE  OF  THE  20 
COSINE  BINS  IN  NRMLZD  IS  SET  UP  ON 
EACH  PASS  THROUGH  DO  14000  I  LOOP. 

I  IS  BIN  SUBSCRIPT  IN  NRMLZD. 

I+l  IS  CURRENT  BAR  IN  ARRAY  LYNPLT. 

DO  14000  I-l.NMUBIN 

hJPLX40  -  40  ♦  NRrt.ZD<  I,  MXAZM.  KTYPE.  KSURF) 

LINTW*  -  NPLX40  /  I  YTOP 

IF  ( <NPLX40-(LINTMP*IYT0P) ) .  GT.  0)  LINTMP-LINTMP  +  l 
LINE  -  42  -  LINTMP 
IF  <LINE.  GT.  1)  GO  TO  5000 

RATIO  »  FLOAT(NRMLZD  (I,  MXAZM,  KTYPE.  KSURF)  )  /  YMAX 
ENCODE  (4,401.  IXCEED(I))  RATIO 
LINE  -  1 
CONTINUE 

JSTART  *  2  <4*<I-1)) 

JSTOP  -  JSTART  +  3 
DO  6000  J-JSTART.  JSTOP 
LYNPLT<J. LINE)  -  IPLUS 
CONTINUE 

IF  (LINE.  GE.  41)  GO  TO  8000 
KSTART  -  LINE  +  1 
DO  7000  K-KSTART.41 

DO  7000  J-JSTART.  JSTOP 
LYNPLT  (J.K)  »  IBLANK 
CONTINUE 

CONTINUE 

IF  (LINE-LAST)  9000.  14000.  I 1000 
C04TINUE 
ICOLUM  =  JSTART 
DO  10000  J-LIhE.  LAST 


LYNPLTdCOLUM,  J)  «  IPLUS 
10000  CONTINUE 

00  TO  13000 
11000  CONTINUE 

ICOLUH  «  JSTART  -  1 
DO  12000  ^MLAST.  LINE 

LYNPLT ( ICOLUH. J)  «  IPLUS 
12000  CONTINUE 

13000  CONTINUE 

LAST  »  LINE 
14000  CONTINUE 

DO  15000  J-LAST.42 

LYNPLT<B1.  J)  -  IPLUS 
15000  CONTINUE 

C  PUT  LEGEND  IN  PLOT.  Cl  MOVE 

C  EXPANDS  LEGEND  LINE  FROM  ARRAY 

C  LEGEND  INTO  PLOT  ARRAY. 

ENCODE  (24.402.  LEGENDd))  SBND(KTYPE) 

CALL  CIMOVE  (LEGEND.  LEGl) 

ENCODE  (24.403. LEGEND(l))  NPRYMS 
CALL  CIMOVE  (LEGEND.  LEG2) 

C  OMIT  OTHER  LEGEND  LINES  FOR  PRIMARY 

C  PARTICLE 

IF  (KTYPE.  GT.  10)  GO  TO  16000 

ENCODE  (24.  404.  LE0END(  1  >  )  AVYELD 
CALL  CIMOVE  (LEGEND.  LEGS) 


a 


I 

g 

g 


» 

I 

» 

I 

J 


ENCODE  (24.  405.  LEOEND(  1  )  )  ERROR 
CALL  CIMOVE  (LEGEND.  LEG4) 

16000  CONTINUE 

C  PRINT  THE  HISTOGRAM 

PRINT  204.  IBLANK.  (LYNPLT ( I.  1  ) .  I-l,  81  ) 

PRINT  204.  (LAflLAX(12-l)>  (LYNOUT(  J,  I ) ,  J«1 . 81 ) .  LABLAX  ( 12-1 ). 
1  (LYNOUT(K.  I).K«82.  324).  I=>1.  10) 

PRINT  204.  LABLAX(l).  (LYNPLT(  I.  42 ) .  I«l.  81 ) 

WRITE  (LUPRT.205)  (POLDEG(  I ).  I-O.  NMUBIN) 

WRITE  (LUPRT.  207)  lYTOP 

RETURN 

END 

BLOCK  DATA  BLOKAA 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


12  MARCH  1986 

BLOKAA  sets  values  for  some  constants  used  in  program  MARPOP. 

Array  dimensions;  values  used  mainly  for  loop  indices. 

MXSURF  s  Maximum  number  of  target  surfaces 
MXTYPE  a  Maximum  number  of  particle  types.  Numbers  1-5 
are  used  for  secondary  particles  (10  were 


permitted  by  MARLOWE  Version 
primaries. 

11).  Type 

1 1 

is 

MXPOLR 

Maximum  number  of  values  for 
angle  bins. 

boundaries 

of 

po  lar 

MXAZM 

•  Maximum  number  of  values  for 
azimuthal  angle  bins. 

boundari es 

of 

COMMON  /MXDMNS/  MXSURF.  MXTYPE.  MXPOLR.  MXAZM 


C 

C 

C 

C 

C 


I/O  Logical  unit  assignments 

LUSPEC  »  Input  parameter  file 
LUPRT  *  Output  print  file 
LUINFO  *  Output  "information"  file 


1 


IS 

s 

k 


1 


o  o  o  o  o  o  o  o  o  o  r>  o  o  o  o  r>  r>  o  o  o  o  r»  o  o  o  o  o  o  o  o  r»  o  o  o  o  ooooooo 


COMMON  AJNXT8/  LU8PEC.  LUPRT,  LUXNFO 


C 

c 


DATA  MX9URF  /2/.  MXTYPE  /llA  MXFOLR  /21/. 
DATA  LU9PEC  /9/.  LUPRT  /4/.  LUXI^O  /9/ 


MXAZM  /21/ 


END 


8UBROUTXNE  CXMOVE  (X80RCE*  XTAROT) 


24  JANUARY  1986 


CXMOVE  EXPANDS  EXQHT  CHARACTERS  FROM  EACH  WORD  OF  ISORCE  INTO 
ARRAY  XTAROT.  ONE  CHARACTER  PER  WORD.  IN  LEFTMOST  BYTE.  WITH  REST 
OF  WORD  FXLLED  WXTH  A8CXX  SPACES  (40  OCTAL.  20  HEXADECIMAL). 


DIMENSION  ISORCEO).  XTAR0T(24> 


DO  1000  X-1.3 
DO  1000  J-1.8 

ITAR0T(8»(l-l)-t'J)  »  X '0020202020202020'  * 

1  AND(SH1FT(180RCE(X}.8*(J-1)>.  X 'FFOOOOOOOOOOOOOO ' ) 

1000  CONTINUE 
RETURN 
END 

SUBROUTINE  FRONT  (KTYPE.  DEEP.  DEPMAX.  LA.  KARMA) 


17  JANUARY  1986 


FRONT  •itractft  and  tabulatas  tha  following  data  for  particlos 
escaping  fro*  tho  front  surfaco:  (a)  ncNsbor  of  atomic  collisions. 

(b)  numbers  of  replacesients  in  replacement  sequences,  and  (c)  maximum 
distance  from  surface.  A  second  entry  point.  FRXNIT.  is  used  to 
initialise  the  arrays  in  which  data  are  tabulated. 


#*»»«  CAU.XNO  PARAMETERS  #«««* 


KTYPE 


( integer. passed)  is  the  particle  type.  MARLOWE  Version  11 
provided  up  to  10  types.  Version  12  provides  only  5  types. 
Therefore,  current  permitted  values  of  KTYPE  are  1-5  and  11. 
for  up  to  9  types  of  secondary  particles  (1-9)  with  a  value 
of  11  indicating  a  primary  particle. 


DEEP 


(real. passed) 
variable  DEEP. 
ALAT(l)  units, 
reached  by  the 


(100  •  (10*z)).  where  d  is 
z  is  the  maximum  depth  (in 


d  ••• 
and 
truncated 
partic  le. 


to  the  next  lowest  0. 


MARLOWE 
MARLOWE 
1  unit) 


DEPMAX 


( real . passed )  Bin  width  (in  MARLOWE  ALAT(l)  units)  for 
maximum  depth  distribution  tables. 


LA 


(  integer,  passed)  m  296*n»  where  m  is  MARLOWE  variable  LA 
(a  particle  type  identifier)  and  n  is  the  number  of  atomic 
collisions  of  the  particle  between  its  deepest  point  in  the 
slab  and  the  escape  (front)  surface. 


KARMA 


»d) 


( integer,  pa 

fc  ♦  69936  ♦  (a  16eb  ♦  ( 16««2)*c  +  (16«*3)«d  +  (16**4)*e) 
where  k  is  MARLOWE  variable  KARMA  and  a.b.c.  d.e  are  the 
numbers  of  replacements  in  each  of  five  replacement 
sequences. 


o  o  o  ooo  ooooooo  ooo  o  n  o  o  n  n  n  o  n  n  n  o  o  n  n  n  n  o  o  o  o  r»  o  o  o  r> 


This  procssting  was  groupsd  togsthsr  and  put  in  a  subroutina  to  B 

ainiaizo  tha  modifications  to  subroutina  CALCN8  in  MARPOP  Varsion  2.0  ^ 

for  TI-A8C. 

Data  ara  tabulatad  in  tha  arrays  MXDEPf  LENBEQ«  and  NCOLIS.  In  £ 

thosa  arrays  tha  diaansion  11  roprasants  typa  of  atom<  passad  to  this 
routina  in  calling  paraaatar  KTYPE.  Only  valuas  1-9  and  11  ara  usad  „ 
by  data  froa  HARLOUE  Varsion  12«  so  naarly  half  of  aach  array  is  not 
usad  (saa  dascription  of  KTYPE  abova.  >  V 


HXDEP  Maximum  dapth  distribution  tabla.  Dapths  ara  maasurad  in 

units  of  MMLOWE  variabla  ALAT(l).  21  bins  ara  providad  for 
aach  particlo  typa.  MXDEP(l<i)  is  count  of  particlos 
raaching  isaxiaua  dapth  lass  than  1  uniti  MX0EP(2ji)  counts 
particlas  raaching  dapth  of  1  unit>  but  lass  than  2  units; 

.  i  MXDEP(20<i)  counts  particlas  raaching  19  units  but 

lass  than  20  units;  MX0EP(21fi)  is  count  of  particlos 
raaching  20  units  or  aora. 

LENSEQ  Langth  of  raplacamant  saquancas  tabla.  19  collisions  in 
saquanco  X  9  saquancos  X  11  particla  typas. 


g 


NCOLIS 


Nufflbar  of  collisions  from  daapast  point  tabla.  First 
runs  1-21  for  1-20  collisions  and  . OT. 20  in  21st  bin. 
Sacond  indos  is  for  11  particla  typa'o. 


indox 

i 


INTEGER  I.  KARMA.  KARTMP.  KC0LI8.  KMXDEP.  KOUNT.  KTYPE.  LA. 
1  LENSEQ.  MXOEP.  NCOLIS.  NSEQ 

REAL  DEEP.  DEPMAX 
DIMENSION  NSEQ (9) 

COMMON  /FRSURF/  MXOEP (21.  11).  LENSEO(  1  9.  9.  1 1 ).  NCOLIS( 21.  1 1 ) 


SAVE 


Noxt  statamant  causas  valuas  of  all  variablas 
on  axit  from  this  subroutina. 


to  ba  savad 

i 


##«*«  Extract  and  bin  maximum  dapth  data  ***** 

KMXDEP  is  first  calculatad  as  numbar  of  DEPMAX  units, 
than  convartad  to  bin  indax. 

1.  OE-8  in  naxt  statamants  pravants  roundoff  arror 

KMXDEP  »  IFIX<  ( (DEEP-AMOD<DEEP.  100.  0)  >/<1000.  OaDEPMAX)  )+l.  OE-8) 
IF  (KMXDEP. OT.  20)  KMXDEP-20 

MXDEP(KMXDEP+l,  KTYPE)  -  1  MXDEP ( KMXDEP+1 .  KTYPE) 


4 


I 


«**««  Extract  and  bin  collision  data  **««* 

KCOLIS  -  LA  /  69936 
IF  (KCOLIS.  OT.  20)  KCOLIS-20 

NCOLIS  (KCOLIS+1.  KTYPE)  -  1  NCOLIS (KCOLIS+1.  KTYPE ) 

*«*««  Extract  and  bin  raplacamant  saquanca  data  »**«* 

KARTMP  -  KARMA  /  296 
KOUNT  -  O 
DO  1000  1-1.9 

NSEQd)  -  M00( KARTMP.  16) 

IF  (NSEQ(I)  .  OT.  0)  KOUNT -KOUNT-t-l 


I 


j 


KARTtf*  »  KARTMP  /  16 
1000  CONTINUE 

00  2000  1-1.5 

IF  (NBEOd).  OT.  0>  LENSEQ(NSEQ(I)/KOUNT,  KTYPE) 
1  1  LEN8EQ(N8EQ(  D.KOUNT.KTYPE) 

2000  CONTINUE 
RETURN 


*««*  *«*«  **«*  *«««  *«*«  *««>«  *«**  *«««  *«««  «««*  «««■«  *«*« 
««*««  »***« 

*****  ENTRY  POINT  TO  INITIALIZE  ARRAYS  ***** 
*****  ***** 

****  ****  ****  ****  ****  *«**  ****  ****  ****  ****  ***«  **** 


ENTRY  FRINIT 


00  3000  1-1/11 
DO  3000  J-1/21 
MXDEP(J, I)  = 
3000  CONTINUE 

DO  4000  1-1/ 11 
DO  4000  J-1/21 
NCOLIS(J/ I) 
4000  CONTINUE 

00  5000  1-1/  11 
DO  5000  J-1/  5 
DO  5000  K-1/ 
LEN8EQ(K/ 

5000  CONTINUE 
RETURN 
END 

SUBROUTINE  INITLZ 


15 

J/  1) 


12  MARCH  1986 


LOGICAL  LOJDEO 
COMMON  /BINS/  POLDEO ( 0 ; 40 > / 
1  AZMOE0<0; 40)/ 

COMMON  /MORBNS/  XMUM1D(20)/ 


QMU(0:40>/  0MU<40)/  NMUBIN. 
AZMRAD(0:  40>/  DAZBIN(40>/  NAZBIN 
PHIMID(40> 


SAVE  values. 


SAVE 


DATA  MXBINS  /40// 


/3.  141592653590/ 


RAOFAC 


****  ****  ****  ****  ****  ****  ****  ***« 


THETA  (POLAR  ANGLE)  BINS 


****  ****  ****  ****  ********  ****  ***« 


Initialize  polar  angle  bin  boundaries:  Non-zero  POLDEG  value 
indicates  input  boundaries  were  in  degrees  and  boundaries  in 
cosine  (QMU)  units  are  calculated  in  DO  3000  I  loop.  If 
POLDEO  values  are  all  zero,  input  boundaries  were  in  mu  units 
and  boundaries  in  degrees  are  calculated  in  DO  4000  I  loop. 


LOJDEO  -  .  FALSE. 


r»  o  o  o  o  o  o  o  o  o  o  o  o  o  o  o  0000000000  000 


DO  1000  I-O.MXBINS 

IF  (POLDEG(  I).  NE.  0.  0)  THEN 
LOJDEO  «  .  TRUE. 

GO  TO  2000 
END  IF 

1000  CONTINUE 
2000  CONTINUE 

IF  (LOJDEO)  THEN 

DO  3000  .I«0.ht1UBIN 

QMU(I)  -  COS(RADFAC»PCLDEO(I ) ) 
3000  CONTINUE 

ELSE 

DO  4000  I-O.NHUBIN 

POLDEGd)  a  ACOS<QMU(I))  /  RADFAC 
4000  CONTINUE 

END  IF 


DO  9000  I-l.NMUBIN 

DMUd)  -  QMUd)  -  QMUd-l) 

XMUMIDd)  a  0.  9  *  (QMUd)  +  QMUd-D) 
9000  CONTINUE 


*  *«««  *«««  *««»  *«■««  »**«  ««*«  «««»  *««-«  » 

*  * 

*  PHI  (AZIMUTHAL  ANCLE)  BINS  * 
»  « 

«  **««  ««*«  *«««  *«««  ««*«  *•»««  ««*«  ««««  « 


Calculate  phi  boundaries  (AZMRAD)  in  radian*  and  width 
(DAZBIN)  and  midpoint  (PHIMID)  of  each  phi  bin. 


DO  6000  laO.NAZBIN 

AZMRADd)  a  RADFAC  *  AZMDEGd) 

6000  CONTINUE 

DO  7000  lal.NAZBIN 

DAZBINd)  a  AZMRADd)  -  AZMRADd-l) 
PHIMIDd)  a  0.  9  ♦  (AZMRADd)  +  AZMRADd-l  >) 
7000  CONTINUE 
RETURN 
END 

SUBROUTINE  INPTA  dl^R*  LUDAT,  KODE) 


24  JANUARY  1986 


INPTA  READS  FIRST  RECORD  WRITTEN  BY  PROGRAM  MARLOW  AND  PRINTS  MARLOW 
RUN  DATE  AND  TIME  IN  HEADING 


*****  CALLING  PARAMETERS  ***** 


ITER  (INTEGER.  PASSED)  FILE  IDENTIFIER 

LUDAT  (INTEGER,  PASSED)  LOGICAL  UNIT  FROM  WHICH  CURRENT  INPUT  FILE 
IS  READ. 

KODE  (INTEGER.  RETURNED)  VALUES  OF  KODE  RETURNED  BY  INPTA: 

1  DATA  READ  OK 

2  End-0 f-file  read,  no  file  FTnn  (nnaLUDAT) 

3  Read  error,  file  FTnn  (nnaLUDAT) 


CHARACTEReS  LDATE.  LTIME.  J0BID(9),  IDENT(19).  OLDTYP.  NEWTYP 


Calculate  width  (DMU)  and  midpoint  (XMUMID)  of  each  theta  bin 


8) 


INTEGER  Z.  ZNREC.  INTEMP.  ITER.  KOOE.  LTYPE.  NTYPE 

ARRAY  INTO  WHICH  RECORD  IS  READ 
DIMENSION  ZNTEMPdSS) 

DIMENSION  LTYPE(IO) 

COMMON  /FSTREC/  ZNREC(128> 

EQUIVALENCE  (NTYPE. ZNREC(7 ) ).  (LTYPE( 1 ) . INREC ( 49) ) 

EQUIVALENCE  (MAXRUN. ZNTEMP (A) ) 

SAVE 

201  FORMAT ('  '.7X. 

1  'PROCESSING  OUTPUT  FROM  PROGRAM  MARLOW  EXECUTION  OF 2  ( IX,  AS)  / 

2  'O'.SAS.  5X.  'MAXRUN-'.  110/'  '.ISAS) 

202  FORMAT< '0'.7X,  '***  UNIT '.  2X,  14,  5X, 

1  'NUMBER  OF  PARTICLE  TYPES  IS',  15,  4X, 

2  'DOES  NOT  MATCH  PREVIOUS:  ',15) 

203  FORMAT  < 'O',  7X,  '***  UNI  T ',  2X,  14,  5X,  'PARTICLE  SYMBOL  NUMBER', 

1  14,  2X,  'IS',  2X,  AS,  4X,  'DOES  NOT  MATCH  PREVIOUS:  ',2X,  AS) 

READ  (LUDAT,END«4000,  ERRS5000)  INTEhF 
CALL  MKCHAR  ( INTEMP (1) , LDATE) 

CALL  MKCHAR  ( INTEMP (3) , LTIME) 

DO  SOO  1-1,5 

CALL  MKCHAR  ( INTEMP  (64-«>I ),  JOBID  ( I ) ) 

SOO  CONTINUE 

DO  900  I-l,  15 

CALL  MKCHAR  <  INTEMP  (69-M  ),  IDENT  ( I ) ) 

900  CONTINUE 

WRITE  (6,201)  LDATE,  LTIME,  JOBID,  MAXRUN,  IDENT 
IF  (ITER.  LE.  1)  GO  TO  3000 

IF  (  INTEMP  ( 7).  EQ.  NTYPE)  GO  TO  1000 

WRITE  (6,202)  LUDAT.  INTEMP (7),  NTYPE 
GO  TO  5000 
1000  CONTINUE 

DO  2000  i-i,hrrYPE 

IF  (  IhfTEMP(4S+I)  .  EQ.  LTYPE(I)  )  GO  TO  2000 
CALL  MKCHAR  (  INTEMP (48-^1 ),  NEWTYP ) 

CALL  MKCHAR  (  IMlEC(48-t>I),  OLDTYP) 

WRITE  (6,  203)  LUDAT,  I,  NEWTYP,  OLDTYP 
GO  TO  5000 
2000  CONTINUE 

3000  CONTINUE 

DO  3500  I-l,  128 

INREC(I)  -  INTEMPd) 

3500  CONTINUE 
KODE  -  1 
RETURN 

4000  CONTINUE 
KODE  -  2 
RETURN 

5000  CONTINUE 
KODE  -  3 
RETURN 
END 

SUBROUTINE  INPT6  (CHI^ODE) 


oooooo  o  o  o  o  r>  r»  o  f>  n  o  o  fi  r»  o  o  o  o  oo  oo  oooooooooo 


12  HARCH  1906 

ZNPTB  raads  parsmatar  racords  10- 

Valuas  of  CHKODE  raturnad  by  INPTB: 

'OK  '  if  data  racords  Mara  raad  aucasafully. 

'EOF  '  if  and-of-fila  raad  (all  racordt  did  not  aiiit). 

'ERROR  '  if  arror  on  attaapt  to  raad  racord. 


CHARACTER»S  FPRIM.  RPRin<  FTARO.  RTARO.  INFO 
CHARACTER  CHK0DE«8.  CHL1NE«81>  CHnUa2 

FolloMing  arrays  provida  for  40  bins#  but  only  20  ara 
alloMad  in  soaia  othar  subprograaiB. 

COMMON  /BINS/  PaL0E0(0: 40) .  QMU(0:40>.  DMU(40)<  NMUBIN* 

1  AZM0E0<0:40).  AZMRAO (0: 40)<  DAZBIN(40).  NAZBIN 

COMMON  /OUTSUT/  FPRIM.  RPRIM.  FTARO.  RTARO.  INFO 
COMMON  /SPECS/  XNORM. 

COMMON  /UNITS/  LUSPEC.  LUPRT.  LUINFO 

SAVE  valuas. 

SAVE 

DATA  MXBINS  /40/ 

100  FORMAT (ABO) 

101  FORMAT  (lOX.  12.  iX.  A2) 

102  FORMAT (10X.7F  10.0) 

103  FORMAT (lOX. FI 0.  0) 

CHKODE  «  'OK  ' 


««««  «»««  «*««  «««»  ««««  ««««  ««#«  «#«« 

*  » 

*  THETA  (POLAR  ANOLE)  BINS  • 

«  « 

*»**  *»*•  ***«  ***«  **«*  ***«  ***«  *•*# 

Raad  numbar  of  mu  (cosina  thata)  bins  to  usa  and  data  typa 
spacifiar  (CHMU).  If  CHMUb'MU'.  raad  bin  boundarias  in 
units  of  mu  and  sat  array  POLDEO  to  zaro.  If  CHMU  is  not 
'MU',  bin  boundarias  ara  in  dagraas.  raad  into  array  POLOEQ. 

Racord  10 

READ  (LUSPEC.  101.  END-3000.  ERR-4000)  NMUBIN.  CHMU 
IF  (CHMU.  EQ.  'MU')  THEN 

♦**  Racord (s)  11 

READ  (LUSPEC.  102.  END-3000.  ERR-4000)  (QMU(  I ) .  I-O.  MIUBIN) 

DO  iOOO  I-O. MXBINS 
POLDEOd)  -  0.0 
1000  CONTINUE 
ELSE 

READ  (LUSPEC.  102.  END-3000.  ERR-4000)  (POLDEG  ( I ) .  I-O.  NfUB  IN) 

END  IF 


»  *«*«  ***«  ♦♦♦*  #*♦*  * 

*  ♦ 

a  PHI  (AZIMUTHAL  ANOLE)  BINS  • 

*  * 

a  *«««  *«*«  ♦*♦♦  *♦♦♦  ***♦  **♦♦  ♦ 


5 

1 

2 
I 
I 
I 
I 
I 
§ 
3 
8 
3 


8 

8 

3 


§ 


1 

s? 


c 

c 

c 


R»ad  ntMibffr  o#  aiiauth^l  anglt  bins  and  bin  boundariss. 


«•«  Racord  12  •** 

READ  (LUBPEC*  101*END«3000.  ERR>4000>  NAZBXN 

*««  Racord(s)  13  »*• 

READ  (LU8PEC«  102.  END-3000.  ERR-4000)  ( AZMDEO(  I ) .  I-O.  KMZBIN) 

»♦♦«»»♦♦♦*#»»>###«#♦#»♦ 
a  * 

•  ROUNDINO  FACTOR  * 


*•*«  ««*«  »«««  ««««  *«««  *«« 


♦**  Racord  14  aa* 


READ  ( LU8P EC.  103.  END-3000.  ERR-4000>  XNORML 

a  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aa 
a  a 

a  OUTPUT  ELIMINATION  SWITCHES  a 
a  a 

a  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aaaa  aa 

Sat  da#aults>  spaci#ving  all  output  is  dasirad. 

FPRIM  -  'FPRIM' 

RPRIM  -  'RPRIM' 

FT  ARC  -  'FT^O' 

RTARO  -  'RTMtO' 

INFO  -  'INFO  ' 

Raad  racord  Fro*  spaci Fication  Fila  and  dacode  data  in  columns 
11-80. 


aaa  Racord  15  aaa 


2000 


3000 


4000 


5000 


READ  (LUSPEC.  100.  ENO-3000.  ERR-4000>  CHLINEd.  BO) 
CHLINE(81:  81)  -  '  ' 

DO  2000  1-11.77 

IF  (CHLIICd:  (Ia4)).  EO.  'FPRIM')  THEN 
FPRIM  -  ' 

ELSE  IF  (CHLINEd:  ( Ia4)  ).  EQ.  'RPRIM')  THEN 
RPRIM  -  ' 

ELSE  IF  (CHLINEd:  <  Ia4)  ).EQ.  'FTARC' )  THEN 
FTARO  -  ' 

ELSE  IF  (CHLINEd:  ( Ia4)  ).  EQ.  'RTARG')  THEN 
RTARO  —  '  ' 

ELSE  IF  (CHLINEd:  (laS)  >.EQ.  'INFO')  THEN 
INFO  -  ' 

END  IF 
CONTINUE 
00  TO  5000 
CONTINUE 
CHKODE  -  'EOF' 

00  TO  5000 

CONTINUE 

CHKODE  -  'ERROR' 

CONTINUE 

RETURN 

END 

SUBROUTINE  MKCHAR  (INWORD.  CHOUT) 


c 

C  23  JANIMRY  1986  11:09 

C 

C  HKCHAR  CONVERTS  THE  VALUE  IN  INUORO.  BYTE>BY-BYTE,  INTO  THE  CHARACTER 
C  VARIABLE  CHOUT. 

C 

CHARACTER*8  CHOUT 
INTEGER  INUORO 
C 

DO  1000  1-1.7 

CHOUT(  1:1)  -  CHAR(  ANO<  X  'OOOOOOOOOOOOOOFF SHIFT  ( INUORO.  8»I )  > ) 
1000  CONTINUE 

CH0UT(8:e)  -  CHAR (AM)<  X 'OOOOOOOOOOOOOOFF'.  INUORO) ) 

RETURN 

END 

SUBROUTINE  OUPTAA 
C 

C  22  MAY  1986 
C 

C  OUPTAA  PRINTS  DISTRIBUTION  TABLES 
C 

CHARACTER*9  FPRIM.  RPRIM.  FTARO.  RTARG 
C  . 

DIMENSION  LABSRF(3.2) 

C 

COMMON  /BINS/  PaLDE0(0: 40) .  QMU(0:  40).  0MU(40).  NMUBIN. 

1  AZMDE0(0:  40).  AZMRAD (0:  40).  DAZBIN(40).  NAZBIN 

COMMON  /FORFAC/  FAC0IV(11.2) 

COMMON  /FSTREC/  ADUMMYO).  MAXRUN.  NTYPE.  BDUMMY(41).  LTYPE(IO). 

1  DDUmY<70) 

COMMON  /MXDf*IS/  MXSURF.  MXTYPE.  MXPOLR.  MXAZM 
COMMON  /MORANS/  XMUMID<20).  PHIMID(40) 

COMMON  /OUTSUT/  FPRIM.  RPRIM.  FTARO.  RTARG 
COMMON  /REAOIN/  SBND(ll).  WIDTH(ll).  DEPBIN.  DEPMAX 
COMMON  /SPECS/  XNORn. 

COMMON  /TABLES/  NENERG  <  100.  11.  2).  NYIELD(21.  10.  2) . 

1  N0EPTH<21.  10.  2).  NRMLZD(21 . 21.  1 1.  2  ) .  FACNRM(  1 1 . 2) . 

2  NAN0LE(21.  21.  11.2).  KSUM<10.2).  KSUMSQ ( 10.  2 ) , 

3  NPRYMS 

COMMON  /UNITS/  LUSPEC.  LUPRT.  LUINFO 
C 

SAVE 

C 

C  Factor  for  converting  degrees  to  radians 

DATA  DEGRAD  /O. 0174932929/ 

DATA  LABSRF  /'  REF',  'LECT'.  'ION  '.  'TRAN'.  'SMIS'.  'SION'/ 

C  MAXIMUM  NUMBER  CF  BINS  IN  'TABLES'  COMMON  BLOCK  ARRAY 

C  NENERG 

DATA  MXMU  /21/.  MXPHI  /21/.  NBIM.M  /lOO/ 

C 

201  FORMAT  < '1 '.  43X.  'PRIMARY  PARTICLE  '.3A4.  '  COEFFICIENTS') 

202  FORMAT  < '1 '.  49X.  3A4.  '  SPUTTERING  OF  '.  A2,  13.  '  ATOMS'/ 

1  IH  .  43X.  'SURFACE  BINDING  ENERGY  (E.  V.  )  '.  G14.  6) 

203  FORMAT < IHO. 92X. 'DISTRIBUTION  OF  YIELD  VALUES'/ 

I  IH  .  8X.  'YIELD'.  8X.  'O '.2019.  '♦') 

204  FORMAT (IH  .  8X.  'FREQUENCY'.  2119 ) 

209  FORMATdHO.  29X.  'DISTRIBUTION  OF  ORIGINAL  DEPTHS  OF  SPUTTERED'. 

1  '  ATOMS  (CHANNEL  WIDTH '.  013.  6.  ')'/ 

2  IH  .8X.  'DEPTH'.  4X.  21 19.  '♦') 

206  FORMAT < IHO. 43X. 'DISTRIBUTION  OF  EJECTED  PARTICLE  DIRECTIONS'/ 


fi  e 


1  '  ',39X,  'POLAR  ANGLE  (ACROSS)'/ 

2  IH  *39X,  'CHANNEL  WIDTH  (DOWN).  AZIMUTHAL  ANCLE  (DECREES)', 

3  Cl  3.  6/ 

4  '  '.  lOX,  'CHAN',  2015,  2X,  'TOTALS'/  ('  ',  lOX,  13,  iX.  2015,  17)  ) 

207  FORMATdHO,  13X,  'PRIMARY  PARTICLES ',  7X,  'MEAN  YIELD',  14X,  'ERROR'/ 

1  IH  ,  18X,  15,  i3X,  FIO.  2,  lOX,  FIO.  2) 

20B  FORMAT  (  'O',  'Total*  (adPhi)  ',2015,  17) 

209  FORMAT ( '1 ', 43X, 'DISTRIBUTION  OF  EJECTED  PARTICLE  DIRECTIONS'/ 

1  'OParticl*  Count/(dMu«dPhi),  normalizad  by  factor  (',C11.4, 

2  '  *  Oil.  4,  '/  ',011.  4,  ')  =  ',011.  4/ 

3  5X.  'Channel ',  IX,  2015/ 

4  9X,  '  dMu  ' ,  3X,  20(  1 X,  F4.  2)  / 

5  2X,  'Phi  ',  5X,  'Mu',  IX,  21(1X,F4.  2)/ 

6  '  Chan  Rad  Deg  I ',  20(4X,  'I '  >,  IX,  'Totals '/ 

7  4X,  F5.  1,  F6.  1,  102X,  '(*dMu) '/ 

8  ('  ',  13,  F5.  1,  F6.  1,2015,  17)) 

210  FORMAT  ('  ',14X,21('I  ')  /  '  POL.  ANO.  DEO ',  21F5.  1  / 

1  '  POL.  ANO.  RAD',  21F5.  1) 

211  FORMAT( 'l'/'0',43X,  'DISTRIBUTION  OF  EJECTED  PARTICLE  DIRECTIONS') 

212  FORMAT  ( 'O',  57X,  'PARTICLE  COUNTS') 

213  FORMAT( 'O',  14X,  'Bin',  2X,20I5  /  15X,  'dMu ',  3X,  20F5.  2  / 

1  16X,  'Mu',  3X,  20F5.  2  / 

2  9X,  'dPhi ',3X,  'Phi',  1X,20(4X,  'I'),4X,  'Row'  / 

3  3X,  'Bin',  1X,2(1X,  '(Deg)'),lX,20(4X,  'I').4X,  'Sums'  // 

4  (4X,  12,  1X,F5.  1,3X,  I3,3X,  2015,  17  / 

5  4(4X,  12,  1X,F5.  1,  3X,  13,  3X. 2015,  17/))) 

214  FORMAT  (8X,  'Column  Sums '.  2X,  201 5,  17  / 

1  'O',  7X,  'Thata  (Deg )',  2X,  2015  / 

2  7X,  'dTheta  ( Deg  ) ',  2X,  20F5.  1 ) 

C 

DELAZM  ■  360.0  /  FLOAT  (NAZBIN) 

AVYELD  »  0  .  0 
ERROR  -  0.  0 
DO  7000  KSURF»1,2 

IF  (((KSURF.  EQ.  D.AND.  (FPRIM.  EQ.  'FPRIM '  ) ).  OR. 

1  ((KSURF.  EQ.  2).  AND.  (RPRIM.  EQ.  'RPRIM')  ))  THEN 

DO  1000  I-1,NBINLM 

IF  (NEhEROd,  11,  KSURF  )  .  NE.  0)  00  TO  2000 
1000  CONTINUE 

END  IF 
GO  TO  3000 
2000  CONTINUE 

C  PRINT  PRIMARY  PARTICLE  TABLES 

PRINT  201,  (LABSRF(  I,KSURF),  I-l,  3) 

KTYPE  -  11 

CALL  OUPTBB  (KTYPE,  KSURF) 

PRINT  206,  DELAZM,  (1,1-1,20), 

1  ((J,  (NANGLE(K,  J,  11,  KSURF),  K-1,  21 )),  J-1,  NAZBIN) 

IF  (NAZBIN.  GT.  1  )  WRITE  (6,208)  (NANGLE(  J,  21 ,  11,  KSURF) ,  J-1,  20) 
IF  (KSURF.  EQ.l)  CALL  OUPTCC  (KTYPE,  DEPM AX) 

CALL  ANOPLT (KTYPE,  KSURF,  AVYELD,  ERROR) 

3000  CONTINUE 

IF  (((KSURF.  EQ.  D.AND.  (FTARO.  EQ.  'FTARO '  ) ) .  OR. 

1  ( (KSURF.  EQ.  2).  AND.  (RTARC.  EQ.  'RTARO')))  THEN 

DO  6000  KTYPE- DhTTYPE 
DO  4000  1-1 ,  NBINLM 

IF  ( NENERO(  I,  KTYPE,  KSURF ).  OT.  0)  00  TO  5000 
4000  CONTINUE 

00  TO  6000 
5000  CONTINUE 


oooooooooooo 


NMUBIN). 


6000 


7000 


PRINT  SECONDARY  PARTICLE  TABLES 
PRINT  202.  (LABSRF(I.KSURF).  I«l,  3).  LTYPE(KTYPE) .  KTYPE. 
i  SBND( KTYPE) 

SUM  «  FLaAT(KSUM(KTYPE.KSURF>) 

PRIMES  «  FLOAT<NPRYMS) 

AVYELO  -  SUM  /  PRIMES 

IF  ( NPRYMS.  LE.  1 )  ERROR  «  O.  O 

IF  <  NPRYMS.  OT.  1 )  ERROR  -  SORT <  (FLOAT  (KSUMSQ(  KTYPE.  KSURF)  )- 
I  (  ( SUM*SUM) /PR  ires )) /FLOAT  ( NPRYMS*  (NPRYMS-1))) 

PRINT  207.  NPRYMS.  AVYELD.  ERROR 
CALL  OUPTBB  (KTYPE.  KSURF) 

PRINT  203.  (1.I-1.20) 

PRINT  204.  (NYIELD(  I.  KTYPE.  KSURF).  1-1.21  ) 

PRINT  209.  OEPBIN.  (I.  1-1.21) 

PRINT  204.  (NDEPTHd.  KTYPE.  KSURF).  1-1.21  ) 

IF  (KSURF.  EO.l)  CALL  OUPTCC  ( KTYPE.  DEPMA X ) 

IF  (NANGLE(MXPOLR.MXAZM.  KTYPE.  KSURF).  OT.  0)  THEN 
WRITE  (LUPRT.  211) 

WRITE  (LUPRT.  212) 

WRITE  (LUPRT. 213)  ( I.  1-1.  NMUB IN) .  (DMU ( J) .  J-1.  NMUB IN) . 

L  ((0.  5*  (QMU  (K)*(1MU(  K-1 )  )).K-1.  NMUBIN). 

1  ( (L. (AZMDEO(L)-AZMDEO(L-l) ). 

}  (NINT(0.  9«(AZM0E0(L)*AZMDE0(L-1))) ). 

)  (NANOLE(M.  L.  KTYPE.  KSURF).  M-1.  NMUBIN). 

)  (NANGLE(MXPaLR.L.  KTYPE. KSURF)) ) . L- 1. NAZBIN ) 

WRITE  (LUPRT.  214) 

L  ((NANOLEd.MXAZM.  KTYPE.  KSURF)).  I- 1 .  htlUB  IN ) . 

2  NANOLE(MXPOLR.  MXAZM.  KTYPE.  KSURF). 

}  ((NINT(0.  9*(P0L0Ee(J)-i>PGLDE0(J-l))  )).  J-1.  NMUBIN). 

)  ( (P0L0E0(K~1)-PQLDE0(K)  ).K-1.  NMUBIN) 

END  IF 

WRITE  (LUPRT.  209)  AVYELO.  XNORML.  FACDIV (KTYPE.  KSURF ) , 
FACNRM(KTYPE.  KSURF).  (I.  I-l.  NMUBIN). 

2  (DMU(J).  J-1.  NMUBIN).  (QMU  (K) .  K-0.  NMUB  IN) . 

2  AZMRAO(O).  AZMOEO(O).  ( (L.  AZMRAD(L).  AZMDEO(L)  > 

)  (NRMLZD(M.L.  KTYPE.  KSURF).  M-1.  NMUBIN). 

)  NRMLZD(MXMU.  L.  KTYPE.  KSURF  ) ).  L-1.  NAZBIN) 

IF  (  NAZB  IN.  OT.  1 )  WRITE  (6.208)  ( NRMLZD(J.  MXPHI.  KTYPE.  KSURF) , 
I  J-1,  NMUBIN) 

WRITE  (6.210)  (POLDEOd).  I-0.IM1UBIN). 

L  (  (DEGRAD«POLDEO(I  )).  I -0.  NMUBIN) 

CALL  ANOPLT( KTYPE,  KSURF.  AVYELD.  ERROR) 

CONTINUE 
END  IF 
CONTINUE 
RETURN 
END 

SUBROUTINE  OUPTBB  (KTYPE.  KSURF) 


4  MARCH  1986 


OUPTBB  PRINTS  ENERGY  DISTRIBUTION  TABLE 


*****  CALLING  PARAMETERS  ***** 


KTYPE  (INTEGER,  PASSED)  PARTICLE  TYPE  (1-10  SECONDARY.  11  PRIMARY) 
KSURF  (INTEGER,  PASSED)  1  -  TARGET  FRONT  SURFACE.  2  -  BACK  SURFACE 


*****  VARIABLE  DECLARATIONS  ***** 


oooo  non  o  oo  r>  n  n  n  n  n  n 


INTEGER  1.  J*  KOUNT.  KSURF.  KTYPE.  NB1N«  NBINLN.  NENERG>  NUN 
REAL  ADUMMY.  WIDTH 
C 

DIHENSIiON  K0UNT(20).  NBIN(20) 

AOUrtlY  IS  NOT  USED  IN  THIS  SUBPROGRAM 
/READIN/  ADUMMY (11).  WIDTH(ll) 

ENERGY  (NENERG).  YIELD  (NYIELD).  DEPTH  (NDEPTH).  AND 
ANGULAR  (NANGLE)  DISTRIBUTION  TABLES. 

KSUM.  KSUMSQ  -  SUM  AND  SUM  OF  SQUARES  OF  YIELD  FOR 

CALCULATING  MEAN  YIELD  AND  STANDARD  DEVIATION. 
NPRYMS  -  TOTAL  NUMBER  OF  PRIMARY  PARTICLES. 

/TABLES/  NENERGdOO.  11.  2).  NYIELD(21.  10.  2) , 

NDEPTH(21>  10.  2).  XANGLE(21 . 21 .  1 1.  2  ) .  FACNRM(11.2). 
NANGLE  (21.  21.  1 1. 2) .  KSUMdO.  2).  KSUMSQ  (10.  2). 
NPRYMS 

/UNITS/  LUSPEC.  LUPRT,  LUINFO 


ENERGY  BIN  DIMENSION  IN  ARRAY  NENERG 
DATA  NBINLM  /I 00/ 

200  FORMAT  ( IHO.  33X.  'EJECTED  PARTICLE  ENERGY  SPECTRUM  (  '.G13.6. 

1  '  E.V.  PER  CHANNEL)'/  IHO.  8X.  'CHANNEL '.  2X.  2015  ) 

201  FORMAT  (IHO.  8X.  'CHANNEL '.  2X .  20(  1 X.  A4 ) ) 

202  FORMAT  (IH  .  8X.  'FREQUENCY'.  2015  ) 

401  FORMAT (14) 

402  FORMAT (13.  '+') 


COMMON 


COMMON 

1 

2 

3 

COMMON 

SAVE 


Print  count*  for  first  20  channsls 

PRINT  200.  WIDTH(KTYPE).  <I.I»1.20) 

PRINT  202.  (NENERGd.KTYPE.KSURF),  I-l.  20> 

Bsyond  channel  20  print  only  non-zero  channels 

Last  channel  is  count  of  "all  greater  than  nex t-to-last" 


NUM  -  0 

DO  2000  I«21.  NBINLM 

IF  ( NENERG (  I.  KTYPE.  KSURF).  LE.  0)  GO  TO  1000 
NUM  -  NUM  +  1 

KOUNT(NUM)  -  hCNERO(  I,  KTYPE.  KSURF) 

IF  (  I.  LT.  NBINLM)  ENCODE(  4.  401.  NB  IN(NUM)  )  I 
IF  (  I.  GE.  NBINLM)  ENCODE(  4.  402.  NB  IN(NUM)  )  I 
1000  CONTINUE 

IF  (  (hWM.  LE.  0).  OR.  (  (NUM.  LT.  20).  AND.  ( I.  LT.  NBINLM)  )  )  GO  TO  2000 
PRINT  201.  (NBIN(J).  J«l.  NUM) 

PRINT  202.  ( KOUNT ( J). J-1. NUM) 

NUM  -  0 

2000  CONTINUE 

RETURN 
END 

SUBROUTINE  OUPTCC  (KTYPE.  DEPMAX) 

C 

C  23  JANUARY  1986 
C 

C  OUPTCC  PRINTS  TABLES  OF  MAXIMUM  PARTICLE  DEPTHS.  NUMBER  OF 
C  COLLISIONS.  AND  NUMBERS  AND  LENGTHS  OF  COLLISION  SEQUENCES  ADDED  TO 
C  OUTPUT  FOR  PROGRAM  MARLOW  VERSION  OF  6  FEBRUARY  1982 


nan  n  o  ooooo 


a 


INTEGER  I.  J.  K.  KTYPE.  L.  LENSEQ.  MXDEP.  NCOL.  NCOLIS.  NROU 
REAL  DEPHAX 

COLUMN  AND  ROW  TOTALS  FOR  COLLISION  SEQUENCE  TABLE 
DIMENSION  NCOLdS).  NR0W<5) 

FRONT-SURFACE-SPUTTERED  PARTICLE  DATA 
MXDEP  «  MAXIMUM  DEPTH  DISTRIBUTION  TABLE 
LENSEQ  •  NUrCER  AND  LENGTH  OF  REPLACEMENT  SEQUENCES  TABLE 
NCOLIS  -  NUmER  OF  COLLISIONS  FROM  DEEPEST  POINT  TABLE. 
FIRST  INDEX  RUNS  1-21  FOR  BIN  VALUES  0-20. 

COMMON  /FRSURF/  MXDEP(  21>  11).  LENSEQ(  1 5.  5.  1 1 ) .  NCOLIS(  21.  1 1 ) 

SAVE 


201  FORMAT ( IHO.  26X.  'DISTRIBUTION  OF  MAXIMUM  DEPTHS  OF  SPUTTERED 

1  'ATOMS  (CHANNEL  WIDTH'.  G13.  4.  ') '/ 

2  IH  .  8X.  'MAX.  DEPTH'.  IX.  2115.  '♦') 

202  FORMAT (IH  .  8X.  'FREQUENCY  '.2115) 

203  FORMAT  ( /IHO.  9X.  'NO.  C^',2SX.  'NUMBER  OF  REPLACEMENTS  IN  SEQUENCE'. 

1  22X. 'NUMBER  OF '/ 

2  IH  .  8X.  'SEQUENCES'.  IX.  1515.  3X.  'PARTICLES'// 

3  5(1H  .  12X.  II.  5X.  1515.  5X.  15/). 

4  IHO.  lOX.  'TOTALS'.  2X.  1515) 

204  FORMAT ( IHO. 39X. 

1  'NUMBER  OF  COLLISIONS  FROM  MAXIMUM  DEPTH  TO  SURFACE'/ 

2  IH  .  BX.  'COLLISIONS'.  4X.  'O'.  2015.  '♦') 


PRINT  201.  DEPMAX.  (I.  1-1.21) 

PRINT  202.  (MXDEP  <  I.  KTYPE).  1-1.21) 

DO  900  1-1 . 5 
NROW(I)  -  0 
900  CONTINUE 

DO  1000  1-1.5 

DO  1000  15 

NROW(I)  -  NROW(I)  LENSEQ(J.  I.KTYPE) 

1000  CONTINUE 

TO  GET  NUMBER  OF  PARTICLES  PER  ROW.  MUST 
DIVIDE  NROW(I)  BY  NUMBER  OF  SEQUENCES  PER 
PARTICLE  IN  ITH  ROW. 

DO  2000  1-2.5 

NROW(I)  «  NROW( I)/I 
2000  CONTINUE 

DO  2500  1-1.15 
NCOL(I)  -  0 
2500  CONTINUE 

DO  3000  1-1.  15 
DO  3000  J-1.  5 

NCOL(I)  -  NCOL(I)  +  LENSEQd.J.  KTYPE) 

3000  CONTINUE 

PRINT  203.  (I.  I-l.  15),  (J,  (LENSEQ(K.J.  KTYPE).  K-1.  15). 

1  NROW(J).  J-1. 5).  (NCDL(L).L-l.  IS) 

PRINT  204.  (I.  I-l.  20) 

PRINT  202,  (NCOLIS(J,KTYPE),  J-1.21  ) 

RETURN 

END 

SUBROUTINE  PRTCLS  ( ISET.  JFYL.  LUDAT.  KODE ) 


C  12  MARCH  1986 
C 

C  On«  MARLOW  'POPDAT'  fil«  (txcvpt  file's  first  record)  is  processed 
C  on  each  coll  to  PRTCLS. 

C  PRTCLS  CALCULATES  NEW  VALUES  FOR  EK  (PARTICLE  ENERGY)  AND  RC0S3 
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("DIRECTION"  COSINE  FROM  NORMAL  TO  TARGET  SURFACE)  AND  COUNTS 
PARTICLES  FOR  ENERGY  (NENERG).  YIELD  (NYIELD),  DEPTH  (NDEPTH), 

AND  POLAR  AND  AZIMUTHAL  ANGLE  (XANGLE)  DISTRIBUTIONS. 

THIS  VERSION  OF  PRTCLS  (FOR  PROGRAM  MARPOP  VERSION  2.  0)  READS  AND 
PROCESSES  RECORDS  FROM  FILE  FT17F001  WRITTEN  BY  PROGRAM  MARLOW 
VERSION  OF  6  FEBRUARY  1982  AND  CALLS  SUBROUTINE  FRONT  TO  TABULATE 
MAXIMUM  DEPTH  DISTRIBUTION  AND  NUMBERS  AND  LENGTHS  OF  DISPLACErCNT 
SEQUENCES. 


•***«  CALLING  PARAMETERS  ****♦ 


ISET  < IntcgeTi  passed > 

JFYL  ( integer/ passed)  Sequence  number  of  current  POPDAT  data  file. 
LUDAT  <  integer/  passed )  Logical  unit  from  uihich  POPDAT  file  is  to  be 
read. 

KODE  (INTEGER.  RETURNED)  VALUES  OF  KODE  RETURNED  BY  PRTCLS; 

1  RECORDS  READ  AND  DATA  PROCESSED  OK. 

2  UNEXPECTED  END-OF-FILE.  NO  DATA. 

3  UNEXPECTED  END-OF-FILE.  PARTIAL  DATA  READ  AND  PROCESSED 

4  READ  ERROR.  NO  DATA  PROCESSED 

5  READ  ERROR.  PARTIAL  DATA  READ  AND  PROCESSED. 

6  MAXIMUM  NUMBER  OF  RECORDS  READ  AND  PROCESSED  NO  EXIT  FLAG 

TEMPORARY  STORAGE  FC«  PARTICLE  YIELD  COUNTS 
DIMENSION  KYELDdOOOO.  10/ 2) 

COMMON  /BINS/  POLDEGCO:  40) .  QMU(0:40).  DMU(40).  NMUBIN. 

1  AZMDEG(0:  40) .  AZMRAD  (0;  40).  DAZB1N(40)/  NAZBIN 

STORAGE  FOR  RECORD  READ  FROM  FILE  FT17F001 
DIMENSION  INREC(8.  16).  XINREC(8,  16) 

ADUMMY  =  ARRAY  NOT  USED  IN  THIS  ROUTINE 
MAXRUN  =  NUMBER  OF  PRIMARY  PARTICLES  IN  MARLOW  RUN 
NTYPE  =  NUMBER  OF  TARGET  PARTICLE  TYPES 
/FSTREC/  ADUMMY  <  5).  MAXRUN.  NTYPE.  BDUMMY(121) 

LIMITS  ON  NRUN  -  NUMBER  OF  PRIMARY  PARTICLES  TO  PROCESS 
FROM  EACH  FT17FXXX  FILE 
/LIMITS/  LIMRUN(IO) 

SBND  -  SURFACE  BINDING  ENERGIES 

WIDTH  =  WIDTH  OF  ENERGY  DISTRIBUTION  BINS  FOR  (1-10) 
SECONDARY  fWD  (11)  PRIMARY  PARTICLES 
DEPBIN  *  WIDTH  OF  BIN  FDR  DEPTH  DISTRIBUTION 
DEPMAX  *  BIN  WIDTH  FOR  MAXIMUM  DEPTH  DISTRIBUTION 
/READIN/  SBND(ll).  WIDTH(ll).  DEPBIN.  DEPMAX 
ENERGY  (NENERG),  YIO-D  (NYIELD).  DEPTH  (NDEPTH).  AND 
ANGULAR  (NANGLE)  DISTRIBUTION  TABLES. 

KSUM.  KSUMSQ  -  SUM  AND  SUM  OF  S(JUARES  OF  YIELD  FOR 

CALCULATING  MEAN  YIELD  AND  STANDARD  DEVIATION. 
NPRYMS  *  TOTAL  NUMBER  OF  PRIMARY  PARTICLES. 

/TABLES/  NENERG  (100.  11.  2).  NYIELD(21.  10.  2) . 

NDEPTH  (21.  10.  2).  X ANGLE (21 . 21 .  1 1.  2  ) .  FACNRM(  1 1 . 2) . 
NANGLE  (21.  21.  1  1.  2) .  KSUM(10.2).  KSUMSQ  ( 10.  2 ) , 
NPRYMS 

/UNITS/  LUSPEC.  LUPRT.  LUINFO 
EQUIVALENCE  ( INREC  (  1.  1  ).  XINREC  ( 1.  1  )  ) 

Save  all  values  on  RETURN  from  this  subroutine 

SAVE 

MAXIMUM  NUMBER  OF  PRIMARY  PARTICLES  PERMITTED  BY  ARRAY 
KYELD 

DATA  MXPTCL  /I 0000/ 


COMMON 


COMMON 


COMMON 


COMMON 

1 

2 

3 

COMMON 


MAXIMUM  NUMBER  OF  BINS  IN  'TABLES'  COMMON  BLOCK  ARRAY 
NENERC 

DATA  NBINLM  /I 00/ 

TWOPI  »  S*PI 

DATA  PI  /3.  141592693590/1  TWOPI  /6.  283189307179/ 


201  FORMAT  ( 'O'.  7X I  '***'.15.  '  PRIMARY  PARTICLES  EXCEEDS  LIMIT  OF',  15. 

1  5X, 'YIELD  FROM  EXCESS  WILL  NOT  BE  COUNTED') 

202  FORMAT  ( 'O'.  7X.  '**♦  NUMBER  OF  AZIMUTHAL  ANCLE  BINS.  '.I5.2X, 

1  'IS  TOO  LARGE.  CHANGING  TO  MAXIMUM  »'.  15) 

203  FORMAT('  BAD  COSINE  IN  ROUTINE  PRTCLS* '.  020.  10,  3X.  'OFF  BY',G20.  10 

1  ',  RESET  TO  1  OR  -1'/ 

2  5X.  'FILE'.  13.  ',  RECORD'.  17.  '.  PARTICLE '.  13,  'DIRECTION  COSINES'. 

3  2(G20.  10,  '.  ').G20.  10) 

204  FORMAT('  MU  OUT-OF-BOUNDS,  SPECIFICATION  SET'.  13. 

1  '.  LOGICAL  UNIT'.  13.  ',  RECORD'.  18,  PARTICLE '.  13,  '.  OFF  BY', 

2  G12.4.  '.  CHANGED  TO'.  F3.0) 

Zero  arrays  before  processing  data  from  first  POPDAT  file 

IF  ( JFYL.  LE.  1 )  THEN 
DO  910  1=1.  2 

DO  910  J=l.  1 1 

DO  910  K=l.  100 

NENERG(K.  J.  I)  =  0 

910  CONTINUE 

DO  920  1=1.  2 

DO  920  J=l.  10 
DO  920  K=l,21 

NYIELD<K.  J,  I)  =0 

920  CONTINUE 

DO  930  1=1,  2 

DO  930  J=l.  10 
DO  930  K-1.21 

NDEPTH<K.  J,  I)  =0 

930  CONTINUE 

DO  960  1=1.  2 

DO  960  J»l.  10 
KSUM(  J,  I)  =  0 
960  CONTINUE 

DO  970  1=1,  2 

DO  970  J=l,  10 

KSUMSQ<J.  I)  =  0 
970  CONTINUE 

DO  985  1=1.  2 

DO  985  J=l,  1 1 
DO  985  K=l,21 
DO  985  L-1.21 

NANGLE<L.  K,  J.  I )  =0 

985  CONTINUE 

NPRYMS  =  0 

Initialize  subroutine  FRONT. 

CALL  FRINIT 
END  IF 

Set  binding  energy  for  primary  particle  and  zero  array 
KVELD  before  processing  data  from  each  POPDAT  file. 


uuuo  oouu  uuao  uuuo  uuouuuu  ouoo  uuuoo 


SBND(ll)  -  0.0 
DO  990  I-1.2 

00  990  ^M1.  10 

DO  990  K>1>  10000 
KYELD(K/  J,  J)  -  0 

990  CONTINUE 


S*t  MAXRUNi  th»  numb»r  of  particlos  to  be  processed  from 
current  POPDAT  file. 

IF  <(LIMRUN(JFYL).  OT.  0  >  .  AND.  (LIMRUN(  JFYL) .  LT.  MAXRUN)  > 

1  MAXRUN  >  LIMRUN(JFYL) 

IF  (MAXRUN.  LE.  MXPTCL)  00  TO  2000 
PRINT  201,  MAXRUN,  MXPTCL 
MAXRUN  -  MXPTCL 
2000  CONTINUE 

AZMBIN  -  TMOPI  /  FLOAT (NAZBIN) 

Each  pass  through  DO  11000  I  loop  processes  one  record 
from  POPDAT  file. 

DO  11000  1-1,1000000 

READ  (LUOAT,  ENO-12000,  ERR-14000)  INREC 

Each  pass  through  DO  10000  J  loop  processes  data  for  one 
partic  le. 

DO  10000  J«l,  16 

INREC(8, J)  (Marlowe  variable  NRUN)  out-of-bounds  is  flag 
that  last  particle  has  been  processed. 

IF  (  (INREC(8,J).  LE.  0).  OR.  (INREC(8.J>  GT.  MAXRUN))  GO  TO  16000 

CALCULATE  EK  CORRECTED  FOR  SURFACE 
BINDING  ENERGY  AND  NEW  Z -DIRECTI ON 
COSINE.  INCREMENT  NENERG  BIN  IF 
PARTICLE  ESCAPES. 

Extract  particle  tgpe  (11-primarg) 

KTYPE  -  MOD(  INREC(7,  J),256) 

IF  (MOD(  INREC(6,  J).2).EQ.  1)  KTYPE-11 

Recalculate  particle  energy  as  (EK»RC0S3e*2?  and  process 
if  particle  escapes  (new  energy  >  SBND) 

TEMP  -  XINRECd,  J)»XIM?EC(4,  J)»XINREC(4,  J)  -  SBNO(KTYPE) 

IF  (TEMP.  LE.  0.  0)  GO  TO  10000 

Calculate  new  direction  cosine  (ZMU)  of  escaped  particle. 
If  calculated  ZI9J  is  invalid  cosine,  adjust  it  and  write 
message  on  file  INFO. 

ZMU  -  SQRT(TEI<«>/(XINREC(  1,  J)-8BN0(KTYPE)  )  ) 

IF  ((ZMU.GT.  1.  0)  .  OR.  (ZMU.  LT.  0.  0)  )  THEN 
IF  (ZMU.  OT.  1.  0)  THEN 
ERRSIZ  -  ZMU  -  1.  0 
ZMU  ■  1.0 
ELSE 
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ERRSZZ  ■  ZMU 
ZMU  -  0.0 
END  IF 

WRITE  (LUIisFO.204)  ISET.  LUDAT.  I,  J.  ERRSIZ.  ZNU 
END  IF 

R*s>t  XINRECdtJ)  to  onorgy  of  oscapod  particle. 

XINREC(1>J>  •  XINREC(1>J}  -  8BND(KTYPE) 

Extract  KARMA  from  INREC(6. j)  and  divide  by  10  to 
determine  surface  from  which  particle  escapes  (KSURF^l 
for  front,  2  for  rear).  KARMA  is  expected  to  be  in 
range  10-29.  Increment  energy  distribution  count. 

KSURF  -  MOD(  INREC(6,  J).69936>/10 

NENBIN  -  1  +  IFIX<XINREC(1.  J)/WIDTH(KTYPE)) 

IF  (NENBIN.  OT.  NBINLM)  hCNBIN-NBINLM 

NENERO<NENBIN,KTYPE,  KSURF)  -  bCNERO(NENB  IN,  KTYPE,  KSURF)  +  1 

If  secondary  particle,  increment  yield  count,  extract 
original  depth  of  particle,  and  increment  depth 
distribution  table. 

IF  (KTYPE.EQ.il)  00  TO  9000 

KYELD(INREC(8,  J),  KTYPE,KSURF)  » 

KYEL0(INREC(8.J),KTYPE,  KSURF)  *  1 
KDEP  «  1  +  IFIX(AM0D(XINREC(9,  J),  100.  0)/DEPBIN) 

IF  (KDEP.  OT.  21)  KDEP>21 

NDEPTH(KDEP,KTYPE, KSURF)  «  NDEPTH( KDEP, KTYPE,  KSURF)  +  1 
CONTINUE 

Calculate  indices  of  mu  and  phi  bins  and  increment 
angular  distribution  table 

DO  6000  KMU«0,NMUBIN 

IF  (ZMU.  OT.  QMU(KMU))  00  TO  6000 
00  TO  7000 
CONTINUE 
KMU  «  NMUBIN  1 
CONTINUE 

FIND  QUADRANT  OF  AZIMUTHAL  ANGLE  AND 
CALCULATE  ANGLE  IN  RADIANS 
COSPHI  «  XINREC(2,  J)/SQRT((XIhFEC(2,  J)*XINREC(2.  J) 

(XINRECO,  J)«XINREC(3,  J)  ) ) 

IF  (  (COSPHI.  OT.  1. 0).  OR.  (COSPHI.  LT.  (-1.  0)  ) )  THEN 
DELTA  -  ABS(COSPHI)  -1.0 

WRITE  (LUIhFO,  203)  COSPHI,  DELTA.  LUDAT.  I.  J. 

XINREC(2,  J),  XIM)EC(3,  J).  XINREC(4.J) 
IF  (COSPHI. OT.  1.0)  THEN 
COSPHI  -1.0 
ELSE 

COSPHI  -  -1.  0 
END  IF 
END  IF 

AZIMTH  -  ACaS(COSPHI) 

IF  (XINRECO,  J).  LT.  0.  0)  AZIMTH  -  TWOPI  -  AZIMTH 
IF  (AZIMTH.  LE.  AZMRAD( NAZ BIN)  )  THEN 
PHI  -  AZIMTH 
ELSE 


PHI  ■  AZIMTH  -  TWOPI 
END  IF 

DO  8000  KPHX-0.  NAZBIN 

IF  (PHI.  OT.  AZMRAD(KPHI)>00  TO  8000 
00  TO  9000 
8000  CONTINUE 

KPHI  •  NAZ6IN  *  1 
9000  CONTINUE 

IF  (  <K«U.  OT.  0).  AND.  (KMU.  LE.  NMUBIN).  AND.  (KPHI.  GT.  0).  AND. 

1  (KPHI.  LE.  NAZBIN) )  MANGLE ( KMU.  KPHI ,  KTYPE,  KSURF  >  »  1  + 

2  NANGLE ( KMU. KPH I . KTYPE.  KSURF ) 

C  TABULATE  MAXIMUM  DEPTH.  NUMBER  AND  LENGTH  OF 

C  DISPLACEMENT  SEQUENCES  AND  NUMBER  OF 

C  COLLISIONS  DATA  FOR  PARTICLES  LEAVING  FRONT 

C  SURFACE 

IF  (KSURF.  EQ.l)  CALL  FRONT  (KTYPE.  XINREC(5.  J).  DEPMAX. 

1  INREC(7.  J).  INREC(6.  J)) 

10000  CONTINUE 

11000  CONTINUE 
KODE  >  6 
GO  TO  17000 
12000  CONTINUE 

IF  (I.  GT.  1 )  GO  TO  13000 
KODE  »  2 
RETURN 
C 

13000  CONTINUE 
KODE  »  3 
GO  TO  17000 
14000  CONTINUE 

IF  (I.  OT.  1 )  GO  TO  15000 
KODE  -  4 
RETURN 
C 

15000  CONTINUE 
KODE  -  5 
GO  TO  17000 
16000  CONTINUE 
KODE  »  1 
17000  CONTINUE 
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ADD  YIELD  DATA  FROM  THIS  FILE  TO 
NYIELD.  KSUM.  AND  KSUMSQ 


DO  18000  KSURF-1.2 

DO  18000  KTYPE-l.NTYPE 
DO  18000  I>1.MAXRUN 
KSUM( KTYPE,  KSURF) 


18000 


NPRYMS  -  NPRYMS  ♦ 

RETURN 

END 

SUBROUTINE  SUMRYZ 


KSUM(  KTYPE,  KSURF)  «  KSUM  (KTYPE.  KSURF)  + 

KYELD(  I.  KTYPE.  KSURF) 
KSUMSQ(KTYPE.  KSURF)  -  KSUMSQ ( KTYPE.  KSURF )  + 

(KYELD(  I.  KTYPE.  KSURF)  ♦  KYELD(  I.  KTYPE.  KSURF) ) 

I BIN  -  KYELD( I. KTYPE,  KSURF)  +  1 
IF  (IBIN.  OT.  21)  IBlN-21 

NYIELD(  IB  IN.  KTYPE.  KSURF)  •  NYIELD  ( IBIN.  KTYPE,  KSURF)  +  1 

CONTINUE 

NPRYMS  *  MAXRUN 


22  MAY  1986 


COHHON  /BINS/  PQLDEC(0: 40> .  QMU(0:40).  DHU(40)>  NNUBIN, 

1  AZMDE0(0:40).  AZMRAD(0: 40)/  DAZBIN(40).  NAZBIN 

COMMON  /FORFAC/  FACD1V(11.2) 

COMMON  /FSTREC/  ADUMMYO).  MAXRUN.  NTYPE/  BDUMMY(121) 

COMMON  /MXOMNS/  MXSURF/  MXTYPE.  MXPGLR/  MXAZM 
COMMON  /SPECS/  XNORM. 

COMMON  /TABLES/  NENEROdOO.  It.  2).  NY1EL0(21/ 10.  2). 

1  N0EPTH(21/  10.2).  NRMLZD(21 . 21.  1 1.  2  ) .  FACNRM(  1 1 . 2) . 

2  NAN0LE(21.  21.  11.2).  KSUMdO.  2).  KSUMSQdO.  2). 

3  NPRYMS 

COMMON  /UNITS/  LUSPEC.  LUPRT.  LUINFO 
C 

SAVE 

C 

C  CALCULATE  ROW  AND  COLUMN  SUMS  FOR 

C  PARTICLE  DIRECTION  TABLE 

DO  7000  KSURF-l.MXSURF 
DO  7000  KTYPE-1. MXTYPE 

IF  <  (KTYPE.  LE.  NTYPE)  .  OR.  (KTYPE.  EQ.  MXTYPE) )  THEN 
DO  2000  I-l.NMUBIN 
NANOSM  «  0 
DO  1900  J-1.  NAZBIN 

NANOSM  »  NANOSM  *  NANOLEd .  J.  KTYPE.  KSURF) 

1500  coivrriNUE 

NANOLEd.  MXAZM. KTYPE.  KSURF)  -  NANOSM 
2000  CONTINUE 

DO  4000  I«l.  NAZBIN 
NANOSM  -  0 
DO  3500  J-l.NMUBIN 

NANOSM  «  NANOSM  <t-  NANOLE(J.  I.  KTYPE.  KSURF) 

3500  CONTINUE 

NANOLE(MXPa.R.  1.  KTYPE.  KSURF)  »  NANOSM 
4000  CONTINUE 

NCROSM  »  0 
DO  6200  I-l.NMUBIN 

NCROSM  -  NCROSM  NANOLEd .  MXAZM.  KTYPE.  KSURF) 

6200  CONTINUE 

NANOLE<MXPaLR. MXAZM.  KTYPE. KSURF)  «  NCROSM 
END  IF 

7000  CONTINUE 
C 

C  Normaliz*  target  values  in  array  NRMLZD  by  factor 

C  (average  yieldelOOOO)  /  total  of  column  with  largest  Mu 

C 

DO  7500  I-l.MXSURF 
DO  7900  J-1. MXTYPE 
DO  7500  K-1.  MXAZM 

DO  7900  L-l.MXPOLR 

NRMLZD(L.  K.  J.  I  )  -  0 

7500  CONTINUE 

DO  9000  KSURF- 1. MX SURF 
DO  9000  KTYPE-1. NTYPE 

IF  ( NANOLE(NMUB IN.  MXAZM.  KTYPE.  KSURF)  .  GT.  0)  THEN 
FACDIV (KTYPE.  KSURF)  -  FLOAT ( 

1  NANGLE(NMUB  IN.  MXAZM.  KTYPE.  KSURF))  /DMU(NMUBIN) 

FACTOR  -  (FLOAT(KSUM{ KTYPE.  KSURF) )  /  FLOAT(NPRYMS ) )  * 

1  XNORML  /  FACDIV(KTYPE.  KSURF) 

FACNRM(KTYPE.  KSURF)  -  FACTOR 
DO  SOOO  I-l.  NAZBIN 
DO  eOOO  J-I.NMUBIN 
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NRHLZDCJ.  I.  KTYPE.K8URF)  •  NINT(FACT0R  *  FL0AT< 
NAN0LE( J. I.  KTYPE. K8URF) ) / (D«U< J> ♦DAZBIN( I > ) ) 

CONTINUE 

DO  8200  I-l.NAZBIN 

NRI1LZD<HXP0LR.  I.KTYPE.  K8URP)  -  NINKFACTOR  *  ( 

FLOAT(NANOLE<HXPaLR. I.KTYPE,  KSURF))  /  DAZBIN(I))) 
CONTINUE 

DO  8400  I-l.NMUBIN 

NRI1LZD<I.HXAZM,KTYPE.KSURF)  -  NINT  (FACTOR  *  < 

FLOAT(NANOLE(I.  MXAZM.KTYPE,  KSURF))  /  DMU<I>>) 
CONTINUE 
END  IF 
CONTINUE 
RETURN 
END 


